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SUMMARY 

Tlie  solution  to  this  problem  has  been  separated  into  two  parts. 

First,  we  find  the  current  distribution  excited  on  the  scatterer  by  an 
external  electromagnetic  wave.  Then,  we  can  calculate  the  radiated  field 
due  to  the  excited  current,  i.e.,  solve  the  Sommerfeld  problem. 

The  earth  is  Idealized  to  be  a homogeneous.  Isotropic  half-space. 

The  scatterers  are  chosen  to  be  perfectly  conducting  coaxial  loops  with 
arbitrary  cross  section  in  the  plane  containing  their  axis.  This  axis 
is  normal  to  the  earth's  surface  and  there  is  rotational  symmetry  so  that 
Fourier  series  expansion  can  be  carried  out  to  simplify  the  problem. 

A variational  approach  has  been  formulated  to  give  the  surface  cur- 
rent distribution  on  the  scatterer.  The  particular  cases  of  scatterers 
consisting  of  one  and  two  cylindrical  loops  (loops  that  are  short  sec- 
tions of  cylindrical  tubes)  have  been  worked  out  explicitly.  Numerical 
computations  still  have  to  be  carried  out  to  see  if  higher-order  terms 
of  the  current  distribution  are  necessary.  Calculations  for  the  current 
distributions  of  loops  with  different  cross  sections  are  desirable  espe- 
cially for  planar  loops,  i.e.,  loops  cut  from  plane  discs. 

Prior  to  obtaining  the  scattered  field,  a comprehensive  review  of 
the  Sommerfeld  problem  was  initiated.  Despite  an  abundance  of  literature 
in  the  area,  the  available  solution  does  not  seem  to  be  conclusive.  This 
review  is  now  in  progress. 

An  experimental  setup  to  check  the  newly  developed  theory  has  been 
devised.  Measurements  of  the  field  scattered  from  a single  loop  above 
the  earth  when  irradiated  by  a 1.5  GHz  electromagnetic  wave  have  been  com- 
pleted. Preparations  are  in  progress  to  carry  out  the  same  experiment 
over  water. 
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EVALUATION 


1.  The  experimental  and  theoretical  work  done  on  this  contract 
shows  how  the  backscattered  energy  from  horizontal  dipoles  and  small 
loops  changes  as  the  characteristics  of  the  earth  change.  (R2B,7.1) 

An  experimental  program  was  carried  out  at  a frequency  of  1.5  GHz. 

The  scattered  field  was  measured  In  elevation  as  the  electrical 
characteristics  of  the  ground  were  changed.  Theoretical  expressions 
for  loops  and  dipoles  near  an  Imperfect  groundplane  were  also  developed. 

2.  This  work  Is  In  support  of  TPO  R2B.  In  order  to  extract  more 
Information  from  the  returned  radar  signal  as  well  as  Improved 
detectability  a better  understanding  of  the  return  from  simple  objects 
must  be  developed.  As  more  low  frequency  radars  (where  the  target 

Is  nearly  resonant)  are  developed  the  characteristics  of  the  earth 
become  more  Important. 

OTHO  E.  KERR 

Project  Engineer 
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1.  INTRODUCTION 


A knowledge  of  the  electromagnetic  field  scattered  by  targets  located 
on  or  above  the  surface  of  the  earth  is  of  primary  Importance  in  radar. 
This  field  is  radiated  by  the  currents  induced  in  the  target  by  the  inci- 
dent plane  wave.  It  follows  that  a first  step  in  evaluating  the  scattered 
field  is  the  determination  of  the  induced  currents.  These  depend  not  only 
on  the  geometrical  and  electrical  characteristics  of  the  scattering  target 
and  the  direction  and  polarization  of  the  incident  wave,  but  also  on  the 
proximity  of  the  earth  with  its  conducting  and  dielectric  properties.  In 
effect,  the  current  in  the  target  corresponds  to  a continuous  distribution 
of  infinitesimal  dipoles,  each  with  appropriate  amplitude  and  relative 
phase.  Thus,  the  total  scattered  field  is  the  vector  sum  of  the  contribu- 
tions by  all  of  the  small  dipoles,  each  with  its  particular  amplitude  and 
phase.  The  field  of  one  of  the  infinitesimal  dipoles  is  the  solution  of 
the  Sommerfeld  problem. 

The  thin  circular  loop  is  the  most  convenient  target  to  analyze 
initially  because  of  its  rotational  symmetry  and  its  usefulness  in  the 
construction  of  composite  targets.  Furthermore,  the  current  distribution 
and  scattered  field  are  well  known  when  the  loop  is  isolated,  i.e.,  suf- 
ficiently far  from  the  earth. 

The  circular  loop  is  also  a simple  target  for  use  in  a systematic 
program  of  measurement.  As  described  in  greater  detail  in  the  Interim 
Scientific  Report,  "Scattering  from  Obstacles  over  the  Earth,"  by  R.  W.  P. 
King,  H.  M.  Lee,  and  L.  C.  Shen  (August  1978)  [1],  one  set  of  measurements 
has  been  completed  of  the  field  scattered  by  horizontal  resonant  and  non- 
resonant circular  loops  and  by  a straight  wire  one-half  wavelength  in 
length  over  a wide  range  of  sizes,  angles  of  incidence,  and  heights  out- 
doors over  the  earth.  That  report  also  contains  a theoretical  treatment 
of  the  scattering  from  horizontal-wire  antennas  over  the  earth. 

An  apparatus  has  been  constructed  to  continue  the  measurements 
systematically  over  a tank  of  water  with  variable  conductivity.  A compre- 
hensive experimental  study  is  about  to  begin.  The  results  should  provide 
an  essential  comparison  with  theory. 


1 

l 


Among  loops  that  possess  rotational  symmetry  around  the  z-axis , which 
Is  normal  to  the  earth'  surface,  the  cylindrical  loop,  which  Is  cut  from 
a cylinder  of  zero  thickness  with  its  height  much  smaller  than  Its  radius 
and  the  Incident  wavelength,  has  the  simplest  integral  equations  for  the 
two  components  of  the  surface  current  excited  by  an  external  field.  The 
two  coupled  Integral  equations  can  be  decoupled,  but  this  Is  of  little 
use.  We  are  faced  with  constants  that  can  assume  very  different  values 
and  affect  the  solution  qualitatively.  The  dielectric  constant  and  con- 
ductivity of  the  earth  vary  from  those  of  dry  sand  to  those  of  sea  water. 

The  loop  can  be  located  close  to  the  surface  or  raised  through  multiples 
of  the  incident  wavelength  until  high  up  in  the  air.  A series  expansion 
of  the  kernel  in  terms  of  the  smallest  scale  factor  for  the  loop,  viz., 
its  height  multiplied  by  the  wave  number  in  the  earth,  may  not  be  a very 
small  number  compared  to  unity.  Hence,  the  expansion  of  the  kernel  in 
terms  of  its  small  parameters  is  a forbidding  task. 

The  Maxwell  equations  for  free  space  are  simple  and  clear.  Rut  in  the 
presence  of  even  the  geometrically  simplest  objects,  they  seldom  permit 
exact  general  solutions  when  subjected  to  the  requisite  boundary  conditions. 
When  applied  to  the  loop,  the  axial  dependence  of  the  current  distribution 
should  be  simple  due  to  the  assumed  smallness  of  the  vertical  dimension  of 
the  cylinder  compared  to  the  wavelength  in  air.  That  is,  we  expect  that 
only  a small  number  of  parameters  whose  magnitudes  are  governed  by  the 
Maxwell  equations  and  the  boundary  conditions  should  adequately  describe 
the  current  distribution.  These  parameters  depend  on  the  size  and  location 
of  the  loop  and  on  the  electrical  properties  of  its  environment. 

A variational  principle  has  been  obtained  that  leads  to  the  coupled 
integral  equations  for  the  components  of  the  current.  Instead  of  attempt- 
ing to  solve  the  integral  equations  directly,  the  variational  integral  is 
used  with  an  approximate  current  distribution  that  is  described  by  the  few 
parameters  we  recognize  as  the  most  important.  This  changes  the  problem 
into  the  relatively  simple  determination  of  several  constants. 

The  variational  approach  is  powerful  not  only  because  it  brings  the 
several  constants  under  our  control,  but  also  because  it  can  be  systemati- 
cally extended  to  multiple  loops,  to  loops  of  different  cross  sections  by 
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the  plane  containing  the  symmetry  axis,  and  to  multiple  loops  each  with 
arbitrary  cross  sections.  This  extension  has  already  been  carried  out. 


2.  FORMULATION  OF  THE  PROBLEM 

‘ 

In  the  cylindrical  coordinates  (p,4,z)  a cylindrical  loop  of  zero 
thickness  is  described  by: 

p-b,  0<h-w<s<h  + v , Ikjwl  « 1 , w/b  « 1 (1) 

The  region  z > 0 is  filled  with  a medium  with  the  complex  permittivity 
ex  and  the  permeability  The  region  z < 0 la  filled  with  a medium 

characterized  by  and  is  assumed  that  e^,  u^,  and  e^,  Uj  are 

constants,  (e ^ + io^/w,  J • 1,  2,  where  o ^ is  the  conductivity.) 

An  external  electromagnetic  wave  E^(p,$,z)  is  incident  on  the 
loop.  It  excites  the  surface  currents  K.($,z)  and  K ($,z)  such  that  the 
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total  tangential  component  of  the  electric  field  vanishes  on  the  loop. 
This  condition  leads  to  the  integral  equation  for  the  surface  current. 
First,  define  the  Fourier  series  expansions  as  follows: 


-E>A/  (b , z) 
▼n 

h +w 

- / b dz 

Wb,z;b,zo) 

G#zn<b»**b»*0) 

* 

«- 

1 

/-s 

N 

O 

• u 

-*(i)(b,z) 

zn 

h-w  0 

G^n(b,r;b,z0) 

G*zn<b»*sb»*0) 

_w. 

where 


<WP’*iV*o>  ' 23?jPp'0  l 1 d> 


(X2  - k2)172 


5“  'Vi<>0> 


Vi^o1  + W>p)Jn+l(iV  - l’2->„(»o).,nOo0) 


V + 


In 


^PPo 


(X2  - k?)1'2 


I - 'l7  “(X2  - k?)2^2(z+  * ) imp 

, J XdX  Rt(X)Jn(Xp)Jn(X»0).  1 


a'  -”Hi 
+ “rr-“ 


l 1 dx  kll(x)J;(x»)J;(xp0).  ' 1 


-(A2-kJ)1/2(z+zQ) 


(7) 


Wp’z;po’zo> 


2a.£l  p 3zQ  / X dA  "(a2  _ ^1/2  Jn(Xp) Jn(XpO) 


[ 


x e 


n2  .2. 1/2.  . 

-(A  -k1)  | z - zQ | 


- Re(A)e 


-(X2-k2)1/2(z  + zn) 


*} 


(8) 


G.*n<p»*ipo*zo> 


l >d> 


1 ^0 


(x2  _ k2j 1/2  Jn(Xp)Jn(Xp0) 


[• 


,.2  . 2v  1/2  i 

• (X  “ K-  ) I 2 • Z 


»'  - Re(X,.-(X2-k2>1/2U+Vj 


C9 


1(011 


1 


(CONTINUED] 


-<»2-kj)1/2l,-g  -(12-kv/2(s+0. 

+ VX>*  J 


[• 


‘“I- 


with 


2 • 

aT3^  l X *“  ~J  . Sjl/i  Jn(Xp)Jn(loO) 
r -(X2  - k2)1/2|.-s  | -<X2-k2)1/2<«+ , n 

L*  - «e(X)«  1 0 J 


R (X)  = 


f2(x2  - k2)1'2  - £-l(x2  - k2)1'2 

e2(>2  - k2)1'2  + c-l(,2  - k2)1'2 


(10) 


(11) 


x2  _ .2.1/2 


R (A)  s 


u,(A‘  - kf)^‘  - uA A"  - kf) 


_ 2 


2 _ .2. 1/2 


..  .2 .1/2  , T2  .2. 1/2 

1*2(1  - k^)  + - k^) 


(12) 


2 2 2 2 

kl  " U Wl^l  * z " 


(13) 


The  branch  cuts  and  integration  path  appear  in  Fig.  1.  Note  that 
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We  know  that  the  functional  V,  viz.. 
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will  give  us  the  Integral  equations  by  means  of  the  calculus  of  variations. 


3.  THE  TEST  FUNCTION 

Because  |k.w|  « 1,  we  expect  K.  (z)  and  K (z)  to  vary  slowly  in  z. 
1 zn 

Furthermore,  we  expect  that  K (h  - w)  - K (h  + w)  - 0,  while 
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We  also  expect  that  the  components  of  the  surface  current  are  well  repre- 
sented by  a few  terms  of  the  following  functions: 
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where  In,  Kn,  s^,  and  t are  parameters  yet  to  be  determined. 

3.1.  The  zero-order  test  function. 

To  lowest  order,  we  assume  that 
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With  this  set  of  current  components,  we  obtain  as: 
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2k, b 


4 1 l lfih  I?'  1 

- TE  l dx  lJ2n(x)  + ln2n(x)1  + [ln  “w"  “ 2 1 

(28) 


Hence,  when  e1  ■ e2,  P2  ■ p2,  we  reproduce  the  result  for  the  theory  of 
the  thin  tubular  loop  antenna  with  width  equal  to  4w  [2], 
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3.2.  The  first-order  test  function. 

We  can  write  the  test  functions  as  follows: 
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The  parameters  *^Qn»  ^nn»  ant*  can  be  obtained  by  solving  the  linear 
equations  in  Table  1 with 
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>>[ 


W(1  - k‘) 


-2h(X2-  kj)1/2 


F2ii(b.w)  - u -^——^0-—  (;)( u° )/ 


J2(*b) 

(X2  - k2) 1/2 


,,2  .2.1/2.  . 

-(X  - kx)  I x - t0[ 


<=„  <»•“)  ■ //  ti  - x2/v2i‘«u . ,2/vV'2  / x dx 

n -w  irV  0 0 (X2  - kj)1/2 

-(X2  - k2)1/2|T-r0| 


^^(b.w.h)  - Jf  ) E^Cb.h+t) 


3 (b,w,h)  - / dr  [1  - t2/w2]1/2  (b , 

n -w  zn 


h + t) 
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Note  that: 

1)  This  variational  principle  provides  the  correct  weight  in  averaging 
the  externally  applied  field  on  the  loop;  and 

2)  The  Integrals  Fq^Oj.w),  F 2 (b#w)t  Gq  (b#w)  can  be  expanded  in 
powers  of  kjW  and  w/b.  After  the  expansion  the  t and  xQ  Integrals  can  be 
easily  Integrated.  This  leaves  the  X Integral  which  can  be  done  numeri- 
cally. 

3.3.  The  lowest-order  test  function  for  two  coupled  loops. 

For  two  cylindrical  loops  described  by: 

pi  ■ bi  » °ihi-wii*iih1  + w1  . IVil  <<  1 <42> 

P2  " b2  . 0 i h2  - w2  i z2  : h2  + «2  • Ik!w2l  K<  1 (43) 

Without  loss  of  generality,  we  may  assume  that  h^  > h2 . The  lowest-order 
test  functions  can  be  chosen  as: 


(*  ~ h) 


2jl72 


(z) 


x[w2  - (z  - h)2]1/2 


(44) 


(45) 


The  variational  Integral  can  be  constructed.  It  leads  to  the  linear  equa- 
tions for  I and  I,  given  in  Table  2,  where 


(A2  - k?)1/2 


^bl*b2 *wi»w2'hl,h2^  " l X dX  pT Re(^)JnUb1)Jn(Xb2) 

« yi.l(A2  - k2)1/2]j0tiv2a2  - k2)1'2,.-''2  * kJ)1/2<hi+V 


(46) 
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^><b1.b2i.1..2,U1.h2)  -lid. 


JOliwl(x2  - kJ)1/2]J0[iw2(A2  - kj)1/2]e  "1 


\(A)j;(Abi)j;(Ab2> 


-(A2-kJ)1/2(h1+h2) 


(47) 


(2) 

F0  /^bl»b2»wl*w2*bl,b2^  " 


,2  “'l  ? "'2 

ii  -i2  .(.2  - 


7,  VXWXV  -(»2-k?)1/2|»1-l.2  + T1-T2 
J A dA  5 o ,/->  e 


, 0 


(A2  - k?)1'2 


(48) 


Note  that  when  ” h2  - wl  + w2* 


2 .2*  1/2, 


(2)  * J„(Xbl)Jn<Xb2>  7 7 

Fon  ^b3.*b2’wi »w2;hi»h2>  - / x dA  "^7^T7r  JoIlwi(A  - 


2 2 1/2  -<x2-‘t5>1/2(h1-h2) 

J0Hw2(XZ  - k2)1/2]e  1 1 2 (49) 


is  a single  Integral. 


4.  DISCUSSION 
(2) 

Except  for  the  constant  ^ (^.^{WW  which  aWears  ln  the 
solution  for  coupled  loops,  each  of  the  coefficients  involved  has  been 
reduced  to  a single  Integral  and  can  be  evaluated  numerically.  Even  Fq2^ 
can  be  reduced  to  a single  Integral  under  the  situations  l^-hj  > w^  + w2 
or  “ b2*  W1  " w2*  Hence»  it  i*  desirable  to  carry  out  the  numerical 
integration  and  make  a comparison  between  the  zero-order  and  first-order 
test  functions  for  the  current  on  a single  loop. 

The  case  of  a single  planar  loop,  i.e.,  one  cut  out  of  a zero  thick- 
ness conducting  disc.  Is  very  Interesting.  To  lowest  order,  the  result 
should  be  the  same  as  thst  of  a cylindrical  loop.  For  higher-order 
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currents  and  for  coupled  loops,  the  answer  is  far  from  obvious. 

We  are  primarily  interested  in  the  scattering  problem.  Now  that 
there  are  test  functions  for  the  currents  it  is  essential  that  we  look 
into  the  effect  of  each  term  of  the  function  on  the  far  field.  A care- 
ful and  detailed  review  of  the  Sommer f eld  problem  is  now  in  progress. 

The  variational  approach  can  be  justified  either  self-conslstently 
by  comparing  results  of  successively  higher-order  test  functions  or  ex- 
perimentally by  comparing  the  calculated  results  with  direct  measurements. 
An  experimental  setup  for  measuring  the  scattered  fields  from  cylindrical 
brass  loops  has  been  completed  {1J.  Unfortunately,  the  C203  generator, 
manufactured  by  the  Applied  Microwave  Division  of  EPSCO,  Inc.,  has  never 
worked  stably  or  delivered  a high  enough  power  output.  It  has  been  deter- 
mined that,  although  the  Siemens  YD1381  tube  used  in  the  generator  may  be 
aging,  the  resonant  cavity  itself  must  first  be  rebuilt  to  correct  a faulty 
design.  The  generator  is  currently  being  repaired. 
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